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Abstract 

The angular coefficients ifk,p( a : 0) of the Fock expansion characterizing the S-state wave function 
of the two-electron atomic system, are calculated in hyper spherical angular coordinates a and 6. To 
solve the problem the Fock recurrence relations separated into the independent individual equations 
associated with definite power j of the nucleus charge Z, are applied. The ’’pure” j-conrponents 
of the angular Fock coefficients, orthogonal to of the hyperspherical harmonics Y^u are found for 
even values of k. To this end, the specific coupling equation is proposed and applied. Effective 
techniques for solving the individual equations with simplest nonseparable and separable right- 
hand sides are proposed. Some mistakes/misprints made earlier in representations of were 
noted and corrected. All j-components of byi and the majority of components and subcomponents 
of ^ 3,0 are calculated and presented for the first time. All calculations were carried out with the 
help of the Wolfram Mathematica. 

PACS numbers: 31.15.-p, 31.15.A-, 31.15.xj, 03.65.Ge 
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I. INTRODUCTION 


The helium isoelectronic sequence presenting a two-electron atomic system contains the 
main features of a many-body system with Coulomb interaction. As such, it can serve as 
a simple basis for testing new quantum theories. The state of a three-body system when 
all the particles are in same space point is known as the triple coalescence point (TCP). 


A long time ago, Bartlett et al. 
be expanded near t 
r 12 . Later Bartlett 
functions 


Q 


showed that the 1 S helium wavefunction, T could not 
re TCP as an analytic series in the interparticle coordinates r i, r 2 and 
2 I and Fock Q| proposed the following expansion containing logarithmic 


[*/ 2 ] 


T(r,a,0) = ^2r k ^2^ p (a,9){lnr) p , 

p =0 


k =0 


where 


a = 2 arctan (r 2 /r 1 ), 9 = arccos \{r\ + r 2 — r5 2 2 )/(2r 1 r 2 )] 


( 1 ) 


( 2 ) 


are the hyperspherical angles, and r = \Jr\ + r\ is the hyperspherical radius. The method 
applied by Fock js| for investigating the 1 S helium wave functions was generalized 0, 5| for 
arbitrary systems of charged particles and for states of any symmetry. The Fock expansion 
(FE) was used for treating the two hydrogen-atom system as the basic one for all subsequent 
calculations in the theory of dispersion forces jjS], The work of Fock was extended by ex¬ 
pansion 1 ^ 2,0 into hyperspherical harmonics (Hff) [7, 8]. The Fock expansion was somewhat 
generalized 9fl to be applicable to any S state, and its first two terms were determined. 
The first numerical solution of the equations for the Fock coefficients was presented in work 


10|] . The most comprehensive investigation on the methods of derivation and calculation of 
the ang ular Fock coefficients was presented in the works of Abbott, Gottschalk and Maslen 
[1 lMl3l| . Methods for simplifying the recurrence relations generated by the Fock expansion 
(TTI) were used ll| to determine the highest power logarithmic terms to sixth order. The 
wavefunction for S states was given to second order in r as a single and double infinite sums 

n n 

111 ]. The results of ref. pL2|] hint at the existence of closed form wavefunction for few-body 
system. The closed form of the helium-like wavefunction including terms up to second order 
in r for 1 S states, and up to forth order for 3 S states was derived in jl3[]. 


In this paper, we build on the work ll|, and therefore we try to adhere to the terminology 
used in this article. We correct some substantial errors/misprints made in the final formulas 
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for ^ 2,0 which is the basis for derivation of the representations for with k > 2. We apply 
a new techniques for calculating an angular Fock coefficients (AFC), and for reducing some 
of them to the form of the one-dimensional series with fast convergence. This technique is 


close to that of used in 


13]. We separate the AFC into the components associated with 


definite power of the nucleus charge, and present all components of " 04,1 and the majority of 
components of for the first time. 

We extensively use all the tools of the Wolfram Mathematica program 
(www.wofram.com/). Its last version, Mathematica 10 will be referred as Mathemat¬ 
ica. To make sure of the correctness of the analytical results, all of them are subjected to 
the numerical verification. 

II. GENERAL APPROACH 


The Schrodinger equation for a system of two electrons, in the field of an infinitely massive 
nucleus, is 

+ = (3) 

where E is the energy, and V is the Coulomb interaction. For a system with nucleus charge 
Z, it is useful to define the dimensionless potential V = Vr, or 

1 


V = 


Z [csc(a/2) + sec(a/2)]. 


( 4 ) 


y/l — sin a cos 9 

The first term of the rhs of Eq. (]4|) presents the electron-electron interaction, whereas the 
second one stands for the electron-nucleus interaction. The Laplacian is 


A = 


1 d 5 d 


-A 2 


where the hyperspherical angular momentum (HAM) operator, projected on S states, reads: 

4 


A 2 = 


sm a 


d 2 d Id d 

2 i tt - sm a — I : — — sm 0 — 
- ' da da sm 6 do do 


( 5 ) 


By substituting the Fock expansion (JT]) into the Schrodinger equation (l3lh one obtains the 
Fock recurrence relation (FRR) 11] 

[A 2 - k(k + 4)] i> k}P = h ktP , 

h k ,p = 2 (k + 2 )(p + l)ijjk,p+i + (p + l)(p + 2) , 0fc,p+2 — 2R'0fc_i iP + 2E^k-2,p- 


3 


(6a) 

(6b) 









Atomic units are used throughout the paper. It is important to note that V’fc.p = 0 for k < 0 


or p > [k/ 2] (see, e.g., 


Ill)- 


We shall solve the FRR (JHJ) to find the angular Fock coefficients, It is well-known 
that the FRR can be solved by expanding the AFC into the hyperspherical harmonics (HH) 
of the form (see, e.g., [ill ]) 

*(t+1) 


Y k i(a,6) = Nki sin* a C*, /2 _)(cos a)P z (cos 9), 


k = 0, 2,4,...; Z = 0,1, 2,..., /c/2 (7) 


where the C”(x) and Pi(z) are Gegenbauer and Legendre polynomials, respectively. The 
normalization constant is 


N kl = 2% 


\2l + l)(k + 2)(k/2-l)\ 
2n 3 (k/2 + l + l)\ ' 


so that 

J Y k i(a,0)Y k 'i'(a, 9)dVL = Skk'Su'i 

where 5 mn is the Kronecker delta, and the appropriate volume element is 


( 8 ) 

(9) 


dQ = 7r 2 sin 2 Q! dasinOdO. a G [0,7r], 9 e [0,7r] (10) 

The HH (J7J) are the eigenfunctions of the operator A 2 , with eigenvalues given by k(k + 4). 
They form a complete set of basis functions in {a, 9}. 

Notice that any general function can be expanded into the hyperspherical harmonics 

T(a,9) = Y / F n iY nl {a,9), ( 11 ) 

nl 

where the expansion coefficients are 

Tni = f 1(«, e)Yra(a, 9)dQ. (12) 

Two important properties of the AFC must be emphasized: 

PI. Any AFC, can be separated into the independent parts (components) 

k—p 

ipkA a ’ 0 ) = ( 13 ) 

j=p 

associated with definite power of Z , according to separation of the rhs (16bj) 

k—p 

= £>£>(<*, 9)Z J (14) 

j=P 
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of the FRR (l6a]i . Hence, each of the FRR (jUJ) can be separated into the individual equations 
(IFRR) for each component: 


[A 2 - k(k + 4)] 6) = h^ p (a, 6). 


O'), 


(15) 


P2. Any component of the angular Fock coefficient, presenting a function of a and 6, must 
be finite at each point of the two-dimensional angular space described by the hyperspherical 
angles a G [0, n\ and 9 G [0,7r]. 

It is well-known that the general solution of the inhomogeneous equation can be expressed 
as the sum of the general solution of the associated homogeneous (complementary) equation, 
and the particular solution of the considered inhomogeneous equation. Note that the linear 
combination of HH, Yki presents the general solution of homogeneous equation associated 
with inhomogeneous Eq. (lT5j) . It must be emphasized that the HH are defined by Eq.(|7]) 
only for even values of k. However, it may be verified that: a) the functions Y^j defined by 
Eq.flZ]), but with odd values of k, also present the solutions of the homogeneous equation 
for Eq. ffl5|l : b) the point a = tt presents a singular one for Yki(a,6) with odd k. It should 
be emphasized that this point presents a pole of the order (/ + 1), and hence, no linear 
combination of Yki with different l (for given odd k) cannot eliminate this singularity. The 
important conclusion is: for odd k the finite particular solution of Eq. fl 15 D presents the 
physical solution we are looking for. The exception is the case of the particular solution, 
which is singular at the point a = tt. This singularity can be removed by subtracting the 
homogeneous solution with equivalent behavior. 

For even k, the solution ^ ctki^ki of the homogeneous equation associated with the IFRR 
(1T51) must be included into the general solution of the IFRR. The coefficients aki for bound 
states are determined by ensuring the wave function is normalisable at r —y oo (see, e.g., 
|l3|). Hence, these coefficients cannot be determined by analysis of the behavior of the wave 
function (jT]) near the triple coalescence point. The exception is the case of k — 2 p, when 


/i 2piP = 0 (see, e.g. jll|). Moreover, it was found |l7j that a 2 o is identically zero (at least) 
for l S state by the required exchange symmetry of the spatial part of the wave function. 
Otherwise, any admixture of Yki preserves the correctness of as the finite solution of 
the IFRR (lisp , making such a solution as the multi-valued one. We present here a single¬ 
valued solutions, which can be produced by orthogonalization of the obtained component 
^kp eac h °f the Tm- The resulting solutions can be called ’’pure” ones, because their HH 
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expansions do not contain Yki for any possible l. 


III. PREVIOUS RESULTS 


In [ll], the angular Fock coefficients were derived for the general Coulomb potential. For 
the case of the helium isoelectronic sequence the following AFC become: 


0i,o(M)) = -Zq+-£, 


•03,1 (a,0) = 

04,2 (a,0) = 


02,1 (a, 9) = -Z 
Z( tt-2) 


7T - 2 

37T 


sin ol cos 9, 


36n 

( 7 r — 2)(57 t — 14) 


180vr 2 


[ 6 Zq sin a cos 9 — £ (6 — 5£ 2 )] , 
Z 2 (l - 2 sin 2 a sin 2 O') , 


where 


r 1 


a 


r 2 


6 = — = cos - , 6 = — = sin - 


a 


S — £1 + € 2 , 


(16) 

(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 


£ = — = Vl — sin a cos 6, 
r 

Using Mathematiaca , we have verified that the AFC presented above satisfy the FRR (l 6 |) . 

The derivation of the AFC, 0fc )O presents the most complicated problem. For k = 2 
this problem was successively solved, first of all, in the works 1 id 13j|. where the S states 
of different symmetry were presented in the natural {ri,r 2 ,ri 2 } coordinates for the general 
Coulomb potential. The expression fo r 0? .n in the, so called, Pluvinage coordinates {£, rj} 
was presented (without derivation) in 15]. In the work jl 6 ] the {0, 77 }-representation from 


was transformed into the {ri,r 2 ,r 12 } coordinates. We have checked both expressions 
for 0 2 ) o, and found that they are, unfortunately, not correct, because they do not satisfy 
neither FRR (J6j) nor the recursion relation expressed in the {£, 77 } coordinates (see, Eqs.(29), 
(25) [3]), whereas the results of 0 do satisfy. Therefore, we present here the closed form 
of the AFC, 02 ,o obtained from the results of jisj], and adapted to the helium isoelectronic 
sequence. The condensed form expressed in the hyperspherical angles {a, 9} reads: 


02 0 = — (1 — 2 E) H- 

12 v ; 6vr 

7 r (y - 4^£) + x/3 


In 


2ny cos 6 hi0 + £) + 7 rx In 
7 — x 


7 + X 


+ i( 2 a 


7T 


+ xa In 


_? 2 (7 + x)_ 

1 + cos 9 


1 — cos 9 


+ 7(2/3 + 7r) + 

+ i xC 




( 22 ) 
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where 


x = cosa, y — sin ci, (23) 

f3 = arcsm(sino;cos 6 l ), 7 = £-\/2 — £ 2 , (24) 

£ = £0 [e i( “-«] + Li 2 [-e- i(ot -^] - Li 2 [—e- i(a+ ^)] - Li 2 [e i( “ +/3) ] . (25) 


Here Li 2 is the dilogarithm function, and i = \/—l. We used the most convenient represen¬ 
tation 


L{4>) = \ [ Li 2 (e 2l< ^) - 1^2(1) - 0(0 - tt)] 


for the Lobachevsky function L( 0 ), which is valid for 0 < 0 < n JL8j, Il9 |. 


In the next sections we shall solve the IFFR (1T51) for the components 0^ of the AFC 
defined by Eq. (fT3|) . We propose special methods for solving Eq. (115]) with different kinds 
of its rhs dehned by Eq. fll 4 [1 and Eq. (16bl) . First, we show that these methods allow us to 
derive the correct expressions for the components of the AFC (I16ll-(fl9ll obtained previously. 
Then, we fix the uncorrect representation for 02,0 obtained in ll|. And at last, we derive 
the components of the AFC 04j and 0 3jO that were not obtained previously. 


IV. TECHNIQUE FOR SOLVING THE IFRR WITH SIMPLEST NONSEPARA- 
BLE RHS 

In this section we discuss the solution of Eq. (1T5]) with the rhs presented by some 
polynomial in variable £ = £(a,0) defined by Eq. (l2Tl) . It follows from Eq. flbb]) . that h i )0 = 
—2U, and hence hf q = — 2/£ presents the simplest example of the rhs mentioned above. We 
shall see that /13 j, /4o, ^ 41 , and many others are the examples of the rhs of that kind. It is 
clear that a physical solution of the corresponding equation (1T5]) reduces to a function $(£) 
of a single variable £. For example, it follows from Eq. (TTfil) that 0 |°q = £/2. It can be shown 
that the result of the direct action of the HAM (J5]) to continues function <f>(£) is 

A 2 4>(£) = (£ 2 - 2 ) $"(£) + ( 26 ) 

Then, Eq. (TT5l) for 0^p(a,0) = *hfc(£) can be rewritten in the form 

(e 2 - 2) m) + + 4 >*»K) = h(o, (27) 
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where h^ p (a,9) = h(£). The general solution of the homogeneous equation associated with 
Eq. (l27|) can be presented in the form 

^fc h) (0 = ClUfc(0 + C 2 Vfc(£), 


where the linearly independent solutions 
Ufc(0 = 


,1/2 

fc+3/2 


h 1 ;:, (£/V2) 


Ql%AUV2) 

v m) = 


(28) 


(29) 


are expressed via the associated Legendre functions, Pj/(x) and Q^{x) of the first and second 
kind, respectively. The particular solution of the inhomogeneous equation (l2Tj) can be found 
by the method of variation of parameters, that yields 

K 




Vk(C)HO 

(£ ,2 -2)W*(£ 




where the Wronskian has a simple form 


w fc (£) = 


+ 2 ) 


(30) 


(31) 


£ 2 ( 2 — £ 2 ) 3 / 2 " 

It was mentioned in Sec. [IT] that the general solution of Eq. (1271) presents a sum of solutions 
([28]) an d (BOD . The coefficients c\ and C 2 in Eq. (BHD mu st be chosen in such a way that 
the final solution becomes the physically acceptable one. It may be verified that Ufc(£) is 
divergent, whereas v fc (£) is finite at the point £ = \[2 (a = tt/ 2,6 = n) for all integer k. 
On the other hand, at the point £ = 0 (a = tt/2,6 = 0) u*(£) is divergent for even values 
of k, whereas v*;(£) is divergent for odd k. This implies the following conclusions. For odd 
values of k, one should put c\ = C 2 = 0 if the particular solution (13H satisfies the finiteness 
condition P2, otherwise the coefficients C\ and/or C 2 must be chosen in such a way that to 
remove the divergence. For even values of k, the additional condition of orthogonality of 
the final solution to Y k i(a,9) enable us to obtain the ’’pure” solutions (see, the end of Sec. 
UT1) . In some complicated cases (see, e.g., Appendix [Cl) the coupling equation (16T1) can be 
applied. 

Eq. flTh]) with all of the right-hand sides mentioned at the beginning of this section, can 
be solved by the method described above. However, for the most of the rhs, one can apply 
the more simple technique described below. 

Substituting <&(£) = H£” into Eq. (l26]h one obtains the following relation 


A 2 BC = BnC [n + 4 — 2(n + 1)£" 2 ] , 









where B is arbitrary constant. Whence, one obtains 


[A 2 - k(k + 4)] BC = BC [(n - k)(n + k + 4) - 2 n[n + l)r 2 ] • (32) 


Given Eq. (l32p . the particular solution (satisfying the finiteness condition) of the correspond¬ 
ing equation 


[A 2 - k(k + 4)] **©= i(0 


(33) 


can be found in the form 


*h ~ 1 

4 >*(£) = ^ 


(34) 


1=0 

where i 0 = 1 for odd n, io — 0 for even n, and ih equals the number of terms in the 
polynomial presenting h(£). The unknown coefficients, Bi can be determined by means of 
substituting (134)1 into Eq.(l33)l. using of Eq . ()32|) . and subsequent equating the coefficients for 
the same powers of £. 

The FRR are solved in order of increasing k and decreasing p. We shall calculate the 
AFC following this rule. Thus, putting ^ 0 ,o = 1, the FRR (IHajl for k — 1, and p = 0 is 
separated into two ones 


(A 2 - 5) V’qo = h i,L U = !) 


,(3) 


(35) 


where hf g = —2/£ (see, the beginning of this section). The use of Eqs. (j32]) - Q3T)) . enables us 
to calculate the component of the AFC, 


= 5 f. 


(36) 


This result is certainly consistent with Eq. dlGjl . Note that all solutions of Eq. dlSp corre¬ 
sponding to the rhs, h(£) (excepting t/>^ which is treated in Sec IVlf All can be calculated 
by the simplified method presented by Eq.( 


V. TECHNIQUE FOR SOLVING THE IFRR WITH SEPARABLE RHS 

This section is devoted to the method of solving Eq. (I15p with the rhs presented by the 
product of functions, each of them depending on only one of the angle variables. 

According to Eq.(j6bJ), the rhs of Eq. (|35]l with j = 1 has a form: 

= 2 [csc(a/2) + sec(a/2)]. (37) 
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A lot of the components of the rhs, presented by Eg. dhbl) . have a form of the product 


h^ p (a,9) = Pi(cos9)(sma) l h.(a), 


(38) 


that for l = 0 reduces to the function of a single variable a, similar to Eq. (l37]h For con¬ 
venience, we have introduced the notation h(a) = h^ p (a). To derive the corresponding 
component of the AFC, we propose the technique described below, 
ft can be shown [13!] that 

d 2 


A 2 Pi (cos 9) f (a) = —4Pi(cos9) 


d 1(1 + 1) 

— + 2 cot a- - 

oa 2 oa sm a 


/(«)• 


(39) 


Hence, the solution of Eq. (fT5]i with i^kp = T)(cos 9)f(a) reduces to Ending the function f(a) 
as a suitable solution of equation 


, d 2 d 41(1 + 1) 

Ittt + 8 cot a o-——- + k(k + 4) 

oa 2 Oa sm" a 


f(a) = — (sina) z h(a). 


(40) 


Putting /(a) = (sina)*g(a), one obtains the following equation 

4g"(ct) + 8(1 + 1) cot a g '(a) + (k — 21)(k + 21 + 4)g(a) = — h(a) 
for the function g(a) = g kp (a). The required solution of Eq. ffIBl) then becomes: 

^p(a,0) = Pi(cos9)(sma) l g(a). 

To solve Eq. (J4TJ) it is convenient to make the change of variable, 


(41) 


(42) 


p = tan(a/2), (43) 

which coincides with definition p — (1 — \x\)/y (given by (All) Q) for 0 < a < 7t/ 2, 
where x and y are defined by Eq. (I23j) . Turning to the variable p, one obtains the following 
differential equation for the function g(p) = g(a), instead of Eq. fllTf . 

(l + p 2 ) g"(p) + 2p 1 [l + p 2 + 1(1 — p 4 )\ g'(p) + (k — 2l)(k + 21 + 4)g(p) = —h(p ), (44) 

where h(p) = h(ci). Using the method of variation of parameters, one obtains the particular 
solution of Eq. (144 [1 in the form 


11 _ f \ f p u ki(p')h(p') j / t ^ r v kl (p')h(p') 7 , /AKS 

g(p) %W J ( i + f f-Y WlW d i> “ w(p) X (i + fPy-w,(f/) dp ' ( ) 


10 













where the linearly independent solutions of the homogeneous equation associated with 
Eq. (l44]) are: 


uM = p- 2l ~\p 2 + 1)S 


\%+l +2 




k + 3 k 1 2 

-,- l + 1 ;- k -P 

2 ’ 2 5 2 ’ ' 


Dtz(p) = (p + l ) 2 


l+Z+2 




P + 3 k 


2 ; 2 +/ + 2 ;/+ 2 ;_pi 


(46a) 

(46b) 


Here 2 -f\ is the Gauss hypergeometric function. It is important to note that: 
1) The corresponding Wronskian 

2/ + 1 fp 2 + l\ a+1 


W t (p) = 


P 


m 


(47) 


is independent of k ; 

2) For the case of k = 21 , the solution u 2 i,i(p) — 1 (follows directly from Eq. 

The lower limits of integration in Eq. (l45p must be chosen in such a way as to remove 
singularities by subtracting the homogeneous solutions with equivalent behavior. This yields 
p c = 1 for even /c, and p c = 00 for odd k. 

One should emphasize that, in fact, functions (|46l) are presented by a rather simple 
elementary functions. For the particular case of the rhs (1371) corresponding to k — l,p — 
0 , l = 0 , j = 1 , p c = 00 , one obtains: 


m = 2(1 + u M = 1 - 3p2 


P 


Application of the form 


Px /1 + pF 

yields for this case: 


vw(p) = 


P* 


W 1 + p 2 


(48) 


',( 1 ) _ 1 + P 


aljW - 


what certainly corresponds to Eq. (fl 6 p . 


= — [sin(a/2) + cos(ck/2)] , 


(49) 


VI. EXPLICATE SOLUTION FOR 

For the case of k — 2, p — 1, the AFC is presented by Eq. (1X7)1 . Derivation of the AFC, 
-0 k,p with k = 2p is a very simple task. Its simplest solution was detailed in Sec.4.3 of the 
paper fll |. 

According to (J 6 ]), the FRR for k = 2, p = 0 is 

(A 2 — 12)02,0 = 802,i — 21401^ + 2 E. (50) 
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Using Eqs.fUl), (fT6|) and (TT7l) . one can present Eq. dSOl) in components 


(A 2 -12 )4 j i = h { 2 % (j — 0,1, 2) 

where 

hfl — 2E — 1, 

/A 1 ) _ u( la ) i jA 16 ) 

n 2,0 — u 2 fl W u 2,0 ) 

h$o = —4(1 + esc a), 


( 51 ) 


(52a) 

(52b) 

(52c) 


(i a ) 2 [sin(a/2) + cos(a/2)] (2 esc a — 3 cos 6 + 3) 
2,0 3\/l — sin a cos 6 

^(ih) csc(a/2) + sec(a/2) 8(7 t — 2) sin a cos 9 
2,0 3\/l — sin a cos 9 37t 

Using the technique presented in SeclY] one obtains : 


1 8 p 

m 20 \ P ) — - + P ~ - r ~:—2 ; 

p 1 + p z 


V20 ( P ) = 

1 + P Z 


Application of the form ([45]) for k = 2,p = 0, l = 0, j = 0, p c = 1 yields 


(53a) 

(53b) 


(54) 


Kp) = 2E - i, 


dD = ^(l-2B). 


Similarly, for k = 2,p = 0, l = 0, j = 2, p c = 1, one obtains: 


Kp) = 


2(1 + P f 

P 


U (2) — - + 

r2,o — o ^ 


1 1 

= —|— sm a. 


3 1 + p‘ 3 2 


(55) 


(56) 


A. Specific solution for U 


(i) 

2,0 


In Sec. IIII1 we presented (without derivation) the closed form of ijj 2j o obtained from the 


results of the work 


13| . This form defined by Eq. (l22|) is very convenient for presentation 


^ 2,0 itself. However, it includes functions (e.g., dilogarithm function with argument in the 
form of the exponential function) which are too complex for further mathematical processing 
required for derivation of the higher-order angular Fock coefficients. For example, according 
to Eq. dfibl) the rhs, h $ t o of the FRR (jfial) for the AFC i/^o, contains the term — 2Vip2, o- To 
obtain the particular solution of the corresponding FRR, one needs to apply the integral 
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formula (j45|) with integrands containing o- It is clear that the reduction of the required 
integration to the analytic form is impossible, whereas the representation of ^ 2,0 hi the 
form of infinite single/double series enables us to provide this integration. Such series were 
presented in [ll]. The problem is that, unfortunately, there are too many errors/misprints 
in the important final formulas, whereas the methods of its obtaining have been described 
rather superficially or not at all been described. 

In this paper, we rederive these formulas, while describing in detail our method that 
differs from that of jll|. 

It is shown in Appendix [A] that the following representation is valid: 


^ 2 }) = — - [sin(a/2) + cos(a/2)] y/l — sin a cos 9 + Xw( a , #), 


where the HH expansion for the function X 20 reads 


X2o(M) = 


D 


nl 


nl 


(n — 2) (n + 6) 


Y nl (a,0), 


(57) 


(58) 


and D n i is defined by Eq. (IA21jh The prime denotes that n = 2 must be omitted from 
summation. Note that (T58|) requires a double summation (over n and /), which possesses a 
slow convergence. 

In [ll] the following single series representation was proposed : 


X2o{a,0) = E Pi (cos 9) (sin a) 1 a 1 (a) 


(59) 


1 =0 


The problem is that the technique, used for deriving the functions cq, is complex and am¬ 
biguous, whereas the final formulas ((A19), (A22) and (A24) from [11]) were presented with 
errors/misprints. We present here the alternative method to derive cq. 

Suppose that a regular function having the unnormalized HH expansion (TTT]h can be 
presented in the form of infinite single series: 


T (a, 9) = ^ Pi(cos9)Qi(a) 


(60) 


1=0 


Multiplication of the right-hand sides of Eq. (JTTj) and Eq.(l60l) by I 2 v,v = (sin a) v Py (cos 9), 
and subsequent integration over the angular space (HH . enables one to obtain the following 
coupling equation: 


A 21,1 — 


(/ + !)! 


a/ttT (Z + 3/2) 


Q i(a) (sin a) l+2 da. 


(61) 
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We used Eg. (|8|) , the orthogonality equation (J9]) for HH’s, and the well-known formula of 
orthogonality for the Legendre polynomials. Putting Qi(a ) = (sina)V;(a), where crja) is 
included into Eq. (EHJ), then using expansion (158|) and representation (IA21D . and simplifying, 
one obtains instead of the latter relation: 


3(7 — 1)(/ + 1)(/ + 3) / (sino:) 2 * +2 cr;(a)cia 

Jo 

0Fr(/ + 3/2) 
2 l+2 l\ 


1 = 


3 r 2 


l -^,- + l,l+-;l + 2,l + 2;l 
2 2 2 


(62) 


Note that the case of / = 1 cannot be used in Eq. (l62j) . because the term with Yji is excluded 
from the HH expansion (l58]h It is clear that for l = 1 the integral in the lhs of Eq. (l62|) just 
equals zero. 

Application of the relations (I40]l- (j44li and (|62l) enables us (details can be found in Ap¬ 
pendix |B]) to derive the following representations which are valid in the range a G [0,7r/2]: 


cr 0 = 


oq = 


12 

1 

247T 


2 y -) a + |x| [1 + 2ln(|x| + 1)] — y 

y / 


(63) 


P 


9 \ i 

-9p — p 3 ) arctan p -- 

P ) P 2 


7T 


~ P 

p (p 2 - 6p + 3) 
72 


- + 16G 


H— In 
6 


1 + p 2 


(64) 


)—i—1 


&l = 


(1 +P Z ) 


2\l—l 


Ip 3 P 2 t P / + 1 

_(/ + !)(/+ 2) ~T + T+l - 1(1- 1) 


+ 


(Z-2)!r(i±l) 


2^! 


l-l 1 + 3 , 3 2 

2 ' 2 ' 2 ’ ^ 


(l > 2) (65) 


r(; + i)r(i + i) 

where G ~ 9159656 is the Catalan’s constant, and x,y and p are defined by Eq. (|23|) and 
Eq.(l48ll. respectively. For tt/ 2 < a < tt, one needs to replace a by 7r — a, x by — x, and 
p by 1/p. One can optionally to set p = (1 — \x\)/y , which is valid for the whole range 
a G [0,7r]. Note that factor 1/4 is missed in the corresponding expression (A24) from 11]. 
Moreover, one should emphasize that representation (A22) from ll| for cri is not correct, 


because it doesn’t satisfy the inhomogeneous differential equation (1B8lh and doesn’t agree 
with definition (l66|) presented below. 


Using series rearrangement (see, 
obtain another representation 

1 


ui. H) of the double summation in Eq. (j58l) . one can 


cp = 


£ 


D 




6 ^ (2m + l — 1)(2 m + l + 3) 

°l, 1 


C. 


(^+i), 


2m 


x), 


( 66 ) 
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which is valid for any l > 0. The convergence of expansion (1661) is very slow, however, it can 
be used to verify the correctness of Eqs. (l63|b([65]h 

It is seen that the component V-4'o is presented by Eqs. (1571b (1551 . On the other hand, using 
definition (1T5D . one can separate the component 02 q out of the AFC -02,o defined by Eq. (122D . 
We have intentionally put up a sign ’’tilde” over 0, because the components obtained by 
these two different methods are not coincident for any point of the angular space {a, 9}. 
It is because both representations mentioned above are not the ’’pure” solutions (see, the 
end of SecHT]) of Eq. (15T]h In other words, the HH expansions of both components include 
some admixture of the homogeneous solutions Y 2 j (l = 0,1) of Eq. (15T|h It is clear that y 2 o 
cannot contain such an admixture by definition (1551) . Thus, using definitions (1571) and (1171) . 
one obtains the coefficient 

C 2 \ = N 21 I Y 2 \ (a, 9) | — - [sin(a/2) + cos(a/2)] y/l — sin a cos 6 1 dQ = , 

for the unnormalized HH, Y 2 i(a,9) = sin a cos 9 in the HH expansion of 020 • H can be 
verihed that C^ = 0. Thus, one obtains the ’’pure” component in the form 0^0^ = 02*0 — 
C 2 i sin a cos It is clear that another way to obtain the ’’pure” component 02 1 ^ is to 
subtract from (a) dehned by Eq. (l64l) . The same method can be used to obtain 
the ’’pure” component Y 2 () > = 02 X o — sin a cos 0 based on the analytic expression (|22l) . 
Numerical integration yields the following value of the coefficient 

C 2 ?=^/r 2 i(«, S )®-( a .^ = 0.315837352, 

whereas C^ = 0. The ’’pure” components 02 1 ^ an d 02'o^ are certainly coincident for any 
angular point under consideration. 

VII. SOLUTIONS FOR AFC WITH k > 2 

For k = 3 and p — 1, the FRR (J6]) reduces to 

[A 2 - 21] 0 3>1 = -2U02,!, (67) 

where V and 02 ,1 are defined by Eq.(j4j) and Eq. (lT7l) . respectively. According to equations 
(fT3]) - (fT5|) . the considered Eq. (l67|) can be presented in components 

[A 2 - 21] 4 ,%} = ft 3 u j, (3 = 1.2) (68) 
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where 


(69a) 

(69b) 


4(1 = 2B(i - ?), 

hf l = —AB cos0[sin(a/2) + cos(ck/2)], 


with constant B = (tt — 2) /37r, and £ defined by Eq. (l2TT) . Using the technique described in 
SeclIVi one easily find the component 


a;= Bt - i). ( to ) 

To solve Eq.flHHD with j = 2, one can apply the technique described in SecJVl Thus, for the 
case of k — 3, p — 1, l — 1, p c = oo, one obtains: 


u( \ no(P + 1 )^/P 2 + 1 t\ 1 + 14p 2 — 35p 4 35 — 14p 2 — p 4 

h(p) = -2 B - x -, «3,i(p) =- „ n - 7 , v 3,i = _ r, - ( 71 ) 


P 


Application of the formula (jj5j) then yields: 


p 3 y/l + P 2 


35y/r+" 


p- 


B( 1 + p) ( 2 ) . Q B |- . f a 

9{P) = TQyTJ 1 ' ' = 9 SU1 a COS 6 = ~2 r V 2 


+ cos 


a 

2/J 


sin a cos 6. (72) 


The components presented by Eq. flTOl) and Eq. (172l) are certainly consistent with AFC -03,1 
presented by Eq. (lT8j) . 


A. The AFC previously undetermined 

Let us consider the FRR for the following two cases: i) k — 3,p — 0, and ii) k = A,p = 1. 


According to (J6J) one obtains: 

(A 2 - 21) ^ 3) o = 1003,1 - 2U^ 2i0 + 2^1,0, (73) 

(A 2 — 32) ip 4 ( i = 2Aip4 t 2 ~~ 2Vip3,i + 2Eip2, i- (74) 

Using Eqs. (TT3|) - (TT5T) and the AFC previously determined, one can present the FRR ([73]) . 
(1771) in components, as follows: 

(A 2 - 21) iP% = h% (j = 0,1, 2, 3) (75) 

(A 2 - 32) V® =/$, (j = 1,2,3) (76) 
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where 


h^ = 2E^-2V 0 ^> = E^ + 


.,( 0 ) 


(o) 


2E-1 
6 £ ’ 


(77) 


= 2 E4]l + 2 V^l - 2 Vo4fi + 10^3 = 

2 . 5(7r-2)(5^ 3 -6n /1-2P\ cr (1 

- 7 X 20 a, 0 + -A- V + -—— — + 2 - - E K 

£ lo 7 r \ 3 / sm a \3 


(78) 


hfl = 2V 1 $l - 2 V q 4 2 ^ + 10^5 = 

4£ 


4c 2 / 

—-X2o(«, 0 ) — 7 ( 2£ + - 

sm a 3 \ 4 


^4,o — 2 Vi'02,o — 


( 2 ) 


3 sin a 
2 


sin a 5 ( 7 r — 2 ) 

H---c sm a cos 6 1 , 


£ 


37T 


3 V sin a 


+ 3 c, 


/ig = 2S^g - 2Wg = [6(1 - 2B) + (12B - 5)f] , 


(79) 

(80) 
(81) 


fcg = 2V,^;; - 2v„^;; + 2^ = ^ x 




(2) 


',( 2 ) 


7T - 2 


2(5tt - 14) 


4 4 

- sin 2 a H— sin 2 aP 2 (cos 0) 
3 3 


57TC 


5 


3 sin ct 


e 5 + 1 - - 


sm a 


ftg = 2nv>g = 


(2) 


7T — 2 
37T 


„ , a\ /a 

2 + tan ( — J + cot y— 


sin a cos 6. 


£- 


■ (82) 
(83) 


Here, Vo and Vi are defined by relation V — Vo — ZVi and Eq.()H), whereas c, £ and y 2 o are 
defined by Eqs. (l20j) . f[2l]) and Eq. (l59l) . respectively. 

It is seen that the rhs and h^\ present functions of only 4, and the corresponding 
components %j)fl and can be derived by simplified method described in SecITVl (see, 
Eqs. ()32]) - (l3ll) ). The results are presented in Table HI The components and obtained 
by technique described in SecjVl are presented in Table Q] as well. 

Extending separation presented in SecO one can obtain the solutions to the majority of 
subcomponents of the remaining components presented by Eqs. f[75]) . f[76j) . In particular, let 
us perform the following additional separation, 


y 3,0 — 3,0 

, (2) _ , (2a) 

y 3,0 — ‘r 3,0 

( 2 ) _ 


+ ^3,0 + ^3,0 + ^3,0 ^ 

(84) 

+ v4 2 0 } + ^3,0 + ^3,o\ 

(85) 

-ff + ^ + 4 2 ?- 

( 86 ) 
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Subcomponents v4!o^ an d are calculated with using Eqs. (T32|) - (134|) . and are presented 
in Table U together with the corresponding rhs. Subcomponents and , obtained 
by the method described in Sec|V], are presented in Table Q] together with the corresponding 
rhs. Note that expression for presented in Tabled] is correct for 0 < a < tt/2, whereas 
for 7 t/2 < a < n one should replace a by n — a. The remaining rhs are: 


(lc) - 2 ~ - e, C = »(«.»). 


h ( c> - 

u 3,0 ~ 


j (26) 5(vr-2) . 

ri\ n = —--sm a cos V, 


''3,0 


37T 


h {2c) - 

u 3,0 — 


4£ 


3 sin a ’ 


h ( 2d) _ 
' i 3,0 — 


4? 


sm a 


X20(«, 9 ), 


(87) 

( 88 ) 


(2c) _ 8(77 — 2) (577 — 14) . 2 


KJ = 
hff = 

4>1 3tt 


457t 2 


a\ /a 

S m|-j+c°s(- 


sin 2 aP 2 (cos 0 ), 
5 




1 - 


£- 


(89a) 

(89b) 


3 sin 01 "* y sin a y ^ £ 

To calculate the subcomponent one can use the particular solution presented in SeclIVl 
Putting k = 3, h(£) = £ 3 in Eq. (j30]h one obtains: 


*3 W K) = § 

*(P)| 


£ (3 - 5£ 2 ) (4 £ 4 - 10 £ 2 + 5) arcsin (£/y/2) 


6 


5v/2^ 

The problem is that < f> 3 ' f4 (£) is singular at the point £ = y/2 (a = tt/2, 6 = n). This 
singularity can be eliminated with the help of the function u 3 (£) (see, Eq. (]29]l ) presenting a 
solution of the associated homogeneous equation, and having the same kind of singularity. 
Thus, given that 

11 r 


Xi = lim < hl p) (£)\/£ — y/2 = -r-rr, 

3 V ;V 80 x 2 3 / 4 ’ 


x 2 = lim i 73 (£) V £ - ^2 = — 7 = 

y/2i r 


one obtains, finally 

z (ic) ( m 25 (tt-2) 
^3,0 (M) = —^— 

25(tt - 2) 
1447T 


■trio---3(f) 

X 2 

£ (3 - 5£ 2 ) (4 £ 4 - 10 £ 2 + 5) arccos (£/ y/2) 

6 + ” 


■37 


(90) 


Using the technique presented in Sec|Vl one obtains for subcomponents with the rhs 
and h,// defined by Eq. ([ 88 j) and Eq. fl89ap . respectively, 


U (2b) - 

X3,0 — 


(n — 2)(1 + p 2 ) -3 / 2 
288vrp 2 


x 


[a — 2p + 14ctp 2 — 35 p 3 ( 7 r — a + 2) — 35p 4 (a + 2) + (14p 5 + p 7 )(vr — a) — 2p 6 ] cos0, (91) 


18 



































(2c) (7T — 2)(57T — 14) . o , . 

ib\ i =--- -sm 2 aP 2 (cos 6) x 

r4 ’ 1 RKdO-rr 2 ’ 


8640vr 2 

{^ ^ (P 2 - !) (3p 8 + 28p 6 + 178p 4 + 28 p 2 + 3) + 6 p (p 8 + 8p 6 + 8p 2 + l)] | , (92) 

where p is defined by Eq. (l43]h Representations (I9TT) and (192|) are correct, as previously, for 
a G [0, 7t/2]. For tt/2 < a < tt one should replace p by 1/p, and a by 7r — a. At first sight 
it may seem that the rhs of Eq. (l9Tj) diverges as a —> 0. However, it may be verified that it 
is wrong, and what is more achieves zero as a —>■ 0. 

For subcomponent corresponding to the rhs (189bl) we use the representation 


= y Pi (cos 0) (sin a) l t/(a), 


(93) 


1=0 


similar to Eq. (|59l) . Function p(p) = t/(cv) can be presented (details can be found in Appendix 
ICl) in the form: 


Mp) = 


it-2 


(19p 5 + 75p 4 - 60p 3 + 30p 2 + 45p - 45) + 


108vr(p 2 + l) 2 115 
^p 5 — 15p 3 + 15p — arctan p + (3p 4 — 10p 2 + 3) 


247 AG , /p 2 +1' 

-In 

757r 7r 


(94) 


h ( p ) = 


7T - 2 


302400ttp 2 (p 2 + 1) 
735 


(l268p 7 - 2505p 6 + 1960p 5 + 32263p 4 + 18900p 3 + 18305p 2 - 735+ 


p 7 + 20p 5 — 90p 3 + 20p + arctan p — 32p 2 (p 2 — 1) In ^ 


(95) 


a ( p ) = 


7T — 2 


14 


41 - 150G + 


3765?r\ 


128 


1 


p‘+ 


141757Tp 4 ( 7T 

5 (672p 9 - 465p 8 + 760p 7 - 6720p 6 - 5880p 5 + 2520p 2 + 315) + 


128 


525 


—-— ( 3p 7 + 28p 5 + 178p 3 + 28p + arctan p — p 4 In ( ^ + ^ 

128 V PJ V 4 


(96) 


r(p) = ( -r 2 y, +1 t 2 x 


1357r2 3 ^ +1 )A;p 2Z+1 


1-3 


p 


21+2 


y plnp n + py P?nP 2n - plo (p 2 + i) 6 arctan p y PlnP^ 


n=0 


n=0 


n=0 


+ A 2 (l)v4i(p), (l > 3) (97) 


where G is the Catalan’s constant, and function n 4 /(p) is defined by Eq. (146bl) for k = 4. The 
A 2 (Z)-factors are presented by Eqs. flC29|) - (lC3ip . whereas the coefficients (i = 1, 2, 3) and 
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A i can be found in Tables munn It is important that A 2 {1) are equal to zero for odd /. Note 
that for / = 1 we have obtained the admixture coefficient T 41 = 0 (see, Appendix IC|) . This 
means that for all odd / in expansion (jU31) . the solution of the form (ITj]) just gives the correct 
results for T\ satisfying the coupling equation (RJD or the equivalent Eq. (|C15jh 

Remind that the general solution of the FRR (jTjj) must contain the addition a^T^, 
where the coefficients au can be determined only by analysis of the asymptotic behavior of 
the wave function (see, the end of Sec lIIIl) . Using Eqs. ([TT]) - ([T2]) . one can detect the presence 
of admixture of the HH Y& in any component U,!m, an d then get rid of such an admixture, 
just as it was done in Sec JVI Al for Therefore, only the ’’pure” components are 

presented by Eqs. (192]) - (jU7|) and Tabled) 

The last subcomponent, we present here, is ■ It is the physical solution of Eq. (1 75 [) 
with the rhs, hf ^ defined by Eq. fl 88 l) . To apply the technique described in SecJV) we again 
use the single sum representation, 

OO 

V4(o C) (M) = ^2 Pi(cos 9)(sina ) 1 (98) 

1=0 

where p is dehned by Eq. (l43lh Function 4>i{p) is obtained (details can be found in Appendix 
iDl) in the form: 

$1 (p) = <t>? (p) + c t v 31 (p), (99) 


where the reduction of the formula (I46bl) for k 


3 and 0 < p < 1 yields: 


V3l{p) = (p 2 + l) 1 1 


{21 — 3)(2/ — 1) 4 2(2/-3) 2 

(21 + 3)(2Z + 5)^ 2/+ 3 P 


( 100 ) 


The particular solution (pf ' 1 of the corresponding differential equation in p (see, Eq. (ID41) ) 
can be presented in the form, 


$\P ) 


2 ~'(p 2 + l) Z ^ 

3(2/ — 3) (2Z — 1) (2Z + 3) (2Z + 5) 


2 fu(p) + 


2 / 2 ; (p) + hi(p) 
2 / + 1 


where 

fu{p) = [9 — 4/(/ + 2)] p + (13 — 4/ 2 ) p 3 , 


f 2 i(p) = [{21 - 3)(2/ - l)p 4 + 2(2/ - 3)(2/ + 5)p 2 + (2/ + 3)(2/ + 5)] arctan(p), 


( 101 ) 

( 102 ) 

(103) 


hi(p) = - [(2/ + 3)(2/ + 5)p 4 + 2(2/ - 3)(2/ + 5)p 2 + (2/ - 3)(2/ - 1)] x 

j + y 2 A] (l, / + 1; / + 2; — p 2 ) . (104) 
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Note that the hypergeometric function presented in Eq. fll()4[) can be expressed through 
elementary functions (see, Eq. (|D 8 |) ). The coefficient q is deffiied by Eqs. (ID17|) - (lD25|) . It 
is clear that formulas (jU9|) - (jlQ4p are valid for 0 < a < 7t/ 2. For 7 t/2 < a < n, one should 
replace p by 1 /p. 

There are only two subcomponents, ^ 30 ^ and V ; 30 > that were not determined in this 
research. The reason is that the right-hand sides and h^ of the corresponding IFRR 
(see, Eqs. (l 8 T]) -([ 88 |) ) include the function X 20 presented by expansion Q59j) , what complicates 
the calculations and dramatically increases the size of the hnal formulas. 


VIII. CONCLUSIONS 


Solutions of the Fock recurrence relations ([ 6 ]) were used to derive the angular coefficients 
'0fc,p(ct, 9) of the Fock expansion (JT|) describing the S'-state wave function of the two-electron 
atomic system. The hyperspherical coordinates with hyperspherical angles (J2]) were applied. 

The FRR were separated into the independent individual equations (TT5]1 associated with 
definite power j (p < j < k — p) of the nucleus charge Z. The appropriate solutions ijjf* of 
Eq. dlSP present the independent components (of the AFC) deffiied by Eqs. (TT3|l - (1T4]) . 

The property of finiteness at the boundary points of the hyperspherical angular space 
was extensively used for deriving each component. 

The ’’pure” components not containing the admixture of the HH, Y^i are found for even 
values of k. 

A few methods for solving the individual FRR were proposed. Simple technique for 
solving the IFRR with simplest nonseparable right-hand side (114)) was described in Sec lIVI 
whereas effective method for solving the IFRR with separable right-hand side of a specific 
but frequent kind was presented in SeclYl 


Some mistakes/misprints made in references [15], 16| for the closed analytic form of 
" 020 ; and - in reference [ll] for the double and single infinite series representations of the 
component 020 > were noted and corrected. 

The coupling equation (l6Tj) was proposed and applied in case of a single series represen¬ 
tation for the component of the AFC. 

Using the techniques mentioned above, all the components of the AFC, 0 41 and the 
majority of components and subcomponents of 0 3O were calculated and presented for the 


21 
















first time in Tables mm and in explicit formulas of Sec JVlll Details of all the calculations 


have been placed in the Appendices lAllDl 

All kinds of calculations (both analytical and numerical) were carried out with the help 
of the program Mathematica (Wolfram.com). 

IX. ACKNOWLEDGMENT 

This work was supported by the Israel Science Foundation grant 954/09 and the ”Pazy 
Foundation”. 

Appendix A 

For ease of comparison, let us present F '20 defined by Eqs. (l5TD . (I52b|) and (155]) . in the 


form 


i> 2 ,o = X 20 (ae,0) + <p(a,0), 


(Al) 


where 


(p(a, 9) = — [sin(a/2) + cos(a/2)] Vl — sin a cos 6. 
3 

By direct action of the HAM operator defined by Eq.(j5]), one obtains: 


(A2) 


(A 2 - 12) <p{a,9) = h [ l 


(A3) 


where the rhs of Eq. (1A3I) is defined by Eq. (l53al) . Hence, 


(A 2 — 12) X 20 («, 9) — /4*o\ 


(A4) 


where the rhs of Eq. (lA4l) is defined by Eq. (l53bl) . 



First of all, we need to solve Eq. (lA4l) by expanding X 20 hi HH. To this end, following 
let us consider the HH expansion of the function 



(A5) 


where Y n i are the unnormalized HH. It follows from (|7|). (1TT1) and dl2]) that 



(A6) 
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According to the works ini and [20J], the following representation holds for v > —2: 


r = 


7 r 


Y\-vj2) 


Zr.fcose) 


1=0 


T(Z — u/2) ( sin a 

r(z +1/2) 




where 


, . . I v l v 1 3 2 

F iA a ) = 2 F M 2 “ 4 ’ 2 “ 4 + 2’ + 2’ Sm a 


(AT) 


(AS) 


presents the Gauss hypergeometric function, and £ is defined by Eq. ([2Tj) . 

Inserting ((7j) (for unnormalized HH) and (1A7D for u — — 1 into flABjl . and using the 
orthogonality of the Legendre polynomials, one obtains 

7T 2 2 1 - 1 r* 


D n i — N nl 


(sin a ) 2l+1 C^+ 2 -i (cos a)f + (a)da. 


21 + 1 


(A9) 


where 


f + (a) = [cos(a/2) + sin(a/2)] 2 F\ ( l; + \,l; + +J + % sin 2 a 


2 4 2 4 2 

It is important to emphasize the following substantial property. It may be verified that 
integral in the rhs of Eq. flA9l) differs from zero only for even values olri/2 — l. 

The following relation was proved to be correct (see, e.g., Eq.(90) jll|L 


/ + («) = 2^1 


l l + l . 2 \ 1. 

-, ——; l + 1; sin a 1 + - sin a 2 F X 


l + l l 2 

~2~’2 + !; Z + 2 ; sin a 


(A10) 


Inserting (lAlOj) into (IA9[) . one obtains: 


D n i — iV 2 z 


2ol -l 


7^2 


2 / + 1 


F\ + ~F 2 ) , 


where 


Zl = 2 f, (1, hti; l + 1 ; y j dx, 


= I y 2 W +% L’,M 2 f, 


'-1 


—+ 1;( + 2;!/ 2 ] circ. 


(All) 


(A 12) 


(A13) 


We made the change of variable, corresponding to notations (1231) . that is x = cosct, y = 
sin ct. Remind the well-known formula of orthogonality for the Gegenbauer polynomials (in 
terms of x, y ): 


1 2 i+i c (i+i) ( \ C (i+i)( x \ dx = n(n + 2l + l)\ 
y o m [X)u n yxjcLx 2 2l + 1 n\(n + l + l)^.) 2 ^ 


(A14) 
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where 5 mn is the Kronecker delta. In order to apply Eq. (1A14D to reduction of integral (IA12D . 
we propose to use expansion (A9) from [ll] with n = —1, v — l + 1, what gives: 

^2 


y 1 2 A, 


where 


l l + l 2 \ rm2 °° 

2’ 2 ’ ' V ) * ^ r 2 (i + m + 3/2) 


2 ,_(/0 v ,(; + 2m+i)r ( m + 1 /2) Gi(;jm)c (j+i )(i)i 


Gi(l,m) = 3 F 2 


l l+l 11 , 3 

2’~ 2 ; 2"”’ i + m+ 2 ;1 


(A15) 


(A16) 


is the generalized hypergeometric function. Inserting (1A15|1 into (1 A1 2[) . and using (1A14I) . 
one easily obtains: 

Pfm -I- 1 /9V/-Urnll 

(A17) 


I 1 = r(m A 1/2 > (, t” )! Gi 


m!r(/ + m + 3/2) 

where m = n/4 — 1/2. 

In order to apply Eq. (1A14I) to reduction of integral (1A13D . we propose to use expansion 
(A6) from ll| with v = l + 1, what gives: 


2 i4 


'+ 1 ,' + Mu ; d^E H4 )" r| ’" +1/ 2) 


7T ^ (Z + m + 1)! 

m —0 




(A18) 


where 


/ 11 3 

G 2 (l , m) = 3 F 2 ( — + 77i +—, — + m + l,l + Tn + — ]l + m + 2,1 + 2m + 2; 1 ) . (A19) 


Inserting (1A19D into (jA13j) . and using (IA14D . one obtains: 

vr(/ + l)(-l) m (2/ + 2m+l)! 


To = 


r G 2 (l,m), 


(A20) 


2 n+l l\m\(l + 2m + 1)(/ + m + 1)! 

where again m = n/4 — 1/2. Substitution of Eqs. (lA17lh (1A20I) and (J8]) into (lAllj) yields 
hnally: 

2 l (/!) 2 (2m)! f 2(1 + 2m + l)r(m + 1/2) (Z + m)\ 


D„i — 


y/Trm\(2l + 2m + 1 )! 


y/ixT(l + m + 3/2) 

(Z + 1)(—1)"T(7 + m + 3/2) 


G\(l, m)+ 
G 2 (l,m) 


(A21) 


2 2m l\(l +m + 1) 

where m = n/4— 1/2 is a non-negative integer. By direct numerical comparing with definition 
(IA6p . it is easy to make sure that relation flA21jl is correct, whereas formula (93) in [ll] is 
not correct. Note that Eq. I1A21 j) is not a single representation for D n i. We have found, at 
least, two another representations. 
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It should be emphasized that representation (1A21D for D n i is correct only for even n/2 — l, 
otherwise, D n i = 0. In particular, D 20 = 0. It follows from (1A21D that D 21 =4 — 8 / 7 r, what 
means that using Eq. (IA5jl and Eq. (l53b|h one can rewrite Eq. (IA4I) in the form 


(A 2 - 12 ) X 2 o(a, 6) = 2 -Y^D nl Y nl {a, 6), (A 22 ) 

nl 

where the prime indicates that n = 2 is omitted from the summation. Taking into account 
that eigenvalues of the A 2 -operator are given by n{n + 4), one obtains for HH expansion, 
Y,ni X n iY n i of the function y 2 o 


(A 2 - 12 ) X 2 o(«, 0) = ^ X n i(n -2)(n + 6 )Y nt (a, 6). (A23) 

nl 

Comparison of (1A22D and (1A23ji yields 

*”<“■*> = !l^' (n-2K: + 6 ) ^( a ' g )' < A24 » 

Remind, that contribution of the terms with n = 2 must be treated separately as the 
solutions of the associated homogeneous equation. 


Appendix B 


To obtain the analytical representations for 07 , it is necessary to substitute expansion (j59l) 
and representation (IA7I1 for v = — 1 into the lhs and rhs of Eci. (IA4lh respectively. Subsequent 
application of the relations (l40]h(l4T| for k = 2 yields 


00 S(tt _S') 

Pi {cos 6 l )(sin a) 1 [4<j/(q!) + 8 (/ + 1) cot a cr'^a) — 4(/ — 1)(/ + 3)<j;(a)] = —--sin a cos 6 


1=0 


a 

esc | — 
2 


sec 


OO / 


sm a 


1=0 


V 2 


37T 


, l 11 3 3 . 2 1 

iFl1 2 + 4’2 + 4 ; + 2 iS1 " 


(Bl) 


Equating the expansion coefficients for the Legendre polynomials (in cos 9) of the same 
order, one obtains the ordinary differential equations that must be solved using the boundary 
conditions P2 and the coupling equation (j62l) . We shall solve these equations using the p 
variable defined by Eq. (l43]h Given that 


sin a = 2p/(p 2 + 1), csc(a/2) + sec(ct/2) = p 1 (p + l)\/p 2 + 1, (B2) 
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the following relationship will be useful for further consideration: 


(1113 3 V \_J (p 2 +1)^, 0< P <1 

2 1 V2 + 4’ 2 + 4’ + 2’ (p* + 1)2 ) \ (p2 + iy H p -2i- K p > 1 

Special cases of l — 0, / = 1 and l > 2 will be considered. 


(B3) 


l = 0 


Equating coefficients for the Legendre polynomials of zero order (/ = 0) in both sides of 
Eq. flBlj) . and turning to the p variable, one can employ Eq . MM for k = 2 and l = 0. This 
yields 


2o"/W i 2 (d 2 + 1) o' 


(p 2 + l) 2 So (p) + 


p 


S* 0 (p) + 12So(p) = -Mp), 


where the substitution g(p) = S 0 (p ) = cr 0 (a) was applied. Using Eqs. 
for the rhs in the range p G [0,1] (0 < a < vr/2): 

(P + 1)(P 2 + 1) 


h(p) = 


3 P 


(B4) 
, one obtains 


(B5) 


For l = 0 one can employ formula with k = 2, p c = 1, that gives the following solution 
of Eq. (IB4I) for 0 < p < 1: 

Sq ^ = P( 2 (P 2 ~ 1) ln l(p 2 + !)/ 2 ] ~ p(3p + 2) - 1} - (p 4 - 6p 2 + 1) arctan p 


12p(p 2 + 1) 


An alternative representation is: 
1 


cr 0 = 


12 


^2 y -^ a + x [1 + 2 ln(x — 1 )] — y — 2 j- . 0<cc< 7 t/2 


(B7) 


where x and y are defined by Eq. (12T1) . It is clear that for 7 t/2 < a < ^ r, one needs to replace 
p by 1/p, x by — x, and a by 7 t — a. Representation (lB7h is simpler than the corresponding 
one (A19) from [ill which meanwhile is correct only for 0 < a < n/2. 

It is easy to verify that <r 0 presented by Eq. (1B7[) corresponds to the ’’pure” component 
-02^ (see, the end of SecJTT]). because the HH expansion coefficient 

J20 oc a 0 (a)Y 2 o(a, 9 )d{} 


equals zero. 
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2 . 1 = 1 


In this case, equating coefficients for P\ (cos 6) = cos 6 in both sides of Eq. dBip . using 
(IB2D . (1B3|) and simplifying, one obtains the equation 


2(2 — p 2 )(l + p 2 ) c ,, N 8(tt - 2) (l + p)(l + p 2 ) 2 


(i+ p 2 ) 2 ^(p)+ “ v " = 

p 671 bp 


0 < p < 1 (B8) 


where S\(p) = a\{o). Formula (1451) with k = 2, / = 1, p c = 1 yields the particular solution 
to Eq. (1B8I) in the form 


^(p) = 


24vr 


9 




2vr' 


—z H-9 p - p arctan p -- — p +2-— 

P A P ) P 2 3 


p (p 2 - 6p + 3) 1 (1 + p 2 

72 6 


(B9) 


To obtain the ’’pure” (single-valued) solution to Eq. (lB8j) . one needs to calculate the coeffi¬ 
cient of the unnormalized HH expansion of the solution (lB9|h Remind that represen¬ 
tation (IB9D is correct only for 0 < p < 1 (0 < a < vr/2). For p > 1 (tt/ 2 < a < tt), one 
should replace p by 1/p in the rhs of Eq. (lB9|h Thus, using Eqs. (TTTjl - (TT2]l for n = 2, 1 = 1 
one obtains: 


J ~21 = I a { f\a)Y 2l {a,'d)dQ = 64TV| l7 r 

14 — 48G + 7r(l + 12 In 2) 
727T : 


2 „2 /•' s i p) (p)p 4 dp r sin0cos 2 ede= 


(1 + P 2 ) 5 " Jo 


(BIO) 


where G is the Catalan’s constant, and N 2 \ is defined by Eq.(j8]). The final result is: G\ = 

Sf(p) - 7~2i ■ 

The same result can be obtained using Eq. (j62]h which for l = 1 reduces to 

ptt/2 


y 4 ai(a)da = 0, 


(BH) 


where cr±(ct) presents the general solution of Eq. (IB8H satisfying the condition P2 of finiteness. 


3. / > 2 

The required function ai (/ > 2) can be certainly obtained by means of application of 
the particular solution (145]) , and subsequent use of Eq. (l62|) . However, in this case we shall 
apply the following simplified procedure. 
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Equating coefficients for P;(cos 6) in both sides of Eq. (IB 11) (for l > 2), using (1B2|) . flB3j) . 
and introducing the function Qi(p) by means of the relationship 

ai(a) = Si(p) = + 1 ) Z_1( 5Kp) ! ( B12 ) 

one obtains the differential equation 


p(p 2 + 1 )Qi(p) + 2 [(/ — 2 )p 2 + / + l] QHp) — 6(/ — l)p Qi{p) — [p + l)(p^ + 1), (B13) 

which is correct only for 0 < p < 1 (0 < a < tt/ 2). Mathematica gives the following general 
solution to the homogeneous equation associated with inhomogeneous equation (IB 131) : 


Q?\l>) = Cl 2*1 


■|(-l;i + |;V)+C2/>- a - 1 


>F 


(B14) 


Assuming that a particular solution has a form a + bp + cp 2 + dp 3 , and substituting the latter 
one into the lhs of Eq. (1B13l) . one easily find this particular solution in the form: 


Q?(p) = S 


Ip 3 


-£ + 


P 


l+l 


(B15) 


2 L(/ + l)(/ + 2) l l + l Z(Z — 1) ' 

It is seen that one should put C 2 = 0 to satisfy the condition P2 of finiteness. Thus, using 
(IB 12 j) . (IB 141) and (IB 151) we presently obtain: 


Slip) = 


i-i -1 


(1 + P 2 ) 


2\l-l 


lp 3 


P + P 


l + l 

_(Z + 1)(Z + 2) l ' l + l /(/- 1) 

2C*! 2 F, 


+ 


l-l 1+3 , 3 , 

+ j ; r 


(B16) 


where the relationship 

'l-l 1 + 3 


2-C1 


; l + o) y — (i + p ) 2 Fi l — i; l + —p 


was used (see, 7.3.1.54 [22]). Now we can apply the coupling equation (1621) to find the 
constant C\ = Ci(l). To this end, first, we need to calculate explicitly the integral 


p 2i+ V / (a)da = /C 1 (/) + C'i/C 2 (/), 


where, using Eq. (IB16l) one obtains: 


ZCi(Z) = 


ZC 2 (Z) = - 


2 l+3 

~1T 

2~ l f 1 


P 


21+2 


lp 3 




P 


/ + 1 


o (p 2 + l)*+ 4 L(/ + !)(/ + 2) Z l + l Z(Z - 1) 

2(+l 


dp, 


y ~" 1 ' 2-P1 




4T + 


(B17) 

(BIS) 

(B19) 









































Mathematica gives for the integral (lB18jh 


MO = W7 


l + l 


31 [ (Z - 1) (2Z + 3) 


2+X 




1 , 7 

— ,1;/ + —; —1 

2 2 


Z 


(Z + 1)(Z + 3) 

To reduce integral (IB 19j) we applied the expansion (T9) from ll|, which for parametrization 
n = l, v = 1/2 becomes: 


(B20) 


2£+l t? 
/ 2-C1 


/-I Z + 3 


■ i + 5'9 2 ) = 5 r2 ( i + ?)E 


2 y 2 

3 F 2 


(—l) m (4m + 1)T (m + 1/2) 

2 / m!(Z + m + 1)!T (Z — m + 3/2) 

m=0 


X 


1 -J-’ ^ - m, Z + m + 2; 1 ) c££ 2) (x). 


Substitution of the latter representation into Eq. (IB19D . and subsequent application of the 
formula of orthogonality for the Gegenbauer polynomials, yields finally: 

2 -‘V?r (,1 + 1 ) 


K 2(0 3r(^)r(iy)' 

Using Eq. (l62li and Eq. (lB17D . one obtains: 

m 


(B21) 


Ci = 


3(Z-l)(Z + l)(Z + 3)/C 2 (Z)’ 


where 


M) = /C 3 (Z) + 1 - 3(Z - 1)(Z + 1)(Z + 3)/Cx(Z), 


ir (j) _ MM + 3 / 2 ) p 

MO - -- 3-c 2 


2 > - + 1,Z + -;Z + 2,Z + 2;1 


(B22) 

(B23) 

(B24) 


The simplest way for reduction of /C(Z) is to calculate numerically the first entries of the 
sequence, that is to calculate /C(Z) for Z = 2, 3,4, 5. Using Mathematica , one obtains the 
required sequence: —5/2,—7/3,—9/4,—11/5. It is seen that/C(Z) = —(2Z + 1)/Z. Simplifying, 
we have, finally: 


Ci = 


(i _ 2 )!r (i±i) 

2r 0 + 2) r (1 +1) 


(B25) 


Appendix C 

Solution of the FRR (f76|) for subcomponent with j = 2d implies that we need to find a 
suitable solution to the equation 

(A 2 -32)y4 2 ? = /4 2 ?, (Cl) 
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where hfj and 
expansion (1A7j) for v - 


are defined by Eq. (189bl) and expansion (1931) . respectively. Inserting 
= —1,1, 3 into the rhs of Eq. (189b|) . one obtains the following expansion 


hff\a,9) = Pi(cos9)(sma) l hi(a), (C2) 

1=0 

for the rhs of Eq. flCl|) . where 


h/fad = 


IT — 2 


3vr(2/ - 1)2* 

5 

(2 1 — 3) sin a 


sin 


a\ fa 

2 ) + C ° S ( 2 


x 


F},3 («) - ( 1 - 


sm a 


Ffi(a) — (2 1 — l)F; i _ 1 (a) 


(C3) 


and where the hypergeometric function F Lu is dehned by Eq. dASD . Thus, according to 
Eq. (|4B with g(a) = t;(a) and k = 4, and expansions (1931) and (1C2I) . one obtains the 
following differential equation 


4t'/(a) + 8(1 + 1) cot a t [(a) + 4(2 — /)(/ + 4)t/(ct) = —lp(a:) (C4) 


for the function t/(a) dehned by Eq. (l93j) . Transformation of Eq. (lC4D to variable p = 
tan(a/2) yields (see, Eq. (l44]) l 

{l + p 2 ) 2 t['( p) + 2p~ l [1 + p 2 + 1(1 -p 4 )] t[(p) + 4(2 - l)(l + 4)r;(p) = -hi(p), (C5) 


where hi(p) = h;(a), and T;(p) = t/(a). Turning to the variable p in Eq. (lC3j) . and simplifying, 
one obtains the following compact representation for the rhs of Eq.(C5): 


x 


, i , = (7T-2j(p+lHp- + lj 

1{P> 3tt( 2/- l)(2/ + 3)2 i+1 


(C6) 


where 0 < p < 1. One of the possible ways to solve Eq. (IC5D is to take into account that its 
rhs (1061) presents some linear combination of the terms of the form: 


K,i(p) = p n (p + l )(p 2 + 1) Z \ 


(C7) 
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with integer n > —2. The particular solution, T n j(p) of Eq. (1C5H with changing /q(p) for 
h n ,i(p) defined by Eq. (1C7|h reads: 


Tn,l(P ) = 


(p 2 + l)' + V+ 2 


2 / + 1 


2 F 1 ( -, 3 - /; - - /; -p 2 


x 


3 F 2 Q,l + 4,Z + ^;l + f,l + ^;-p 2 )p | 3 F 2 (|, l + 4, l + l + |, l + -p 2 ) 


2/ + n + 4 


2/ + n + 3 


— 2 -E 1 ( l + 4; l + —p 2 


x 


3 F 2 (f, 3 - Z, | - Z, ^;-p 2 )p 3 F 2 & 3 - I, f + 1; | - Z, f + 2; -p 2 ) 


n + 3 


n + 2 


(C8) 


This solution was derived by means of application of Eq. (jj5j) with changing definite integra¬ 
tion for indefinite integration over p (that is putting the antiderivatives for the lower limit 
equal zero). 

Here, we present the particular solutions, r/ p ' ) of Eq. (IC5D obtained by the exact formula 
(145|) with k = 4. We consider the cases of l = 0,1, 2 (for l max = k/ 2), and l > 3, separately. 

Application of Eqs. (l45|b(l47|) with h(p) = hi(p ) defined by Eq. (IC6|h and p c = 1 yields for 

l = 0 , 1 , 2 : 


Tq P \p) = 


7T — 2 


108vr (p 2 + l) z 


(1 ~P) 
15 


(l9p 4 + 142p 3 + 82p 2 - 48p - 3) - 


p° — lSp 13 + 15p-J arctan p + (3p 4 — 10p 2 + 3) In 


p 2 + 1 


(C9) 


t[ p \ P ) = 


7T - 2 


302400vrp 2 (p 2 + 1) 


[l268p 7 - 2505p 6 + 1960p 5 + 32263p 4 + 18900p 3 + 18305p 2 


735 ( p 7 + 20p 5 - 90p 3 + 20p + arctanp - 23520p 2 (p 2 - 1) In ( p + 1 


735+ 

(CIO) 


A f\p) = 


7T — 2 


362880vrp 4 


[672p 9 - 465p 8 + 760p 7 - 6720p 6 - 5880p 5 + 8798p 4 + 2520p 2 + 315+ 


105(p 2 - 1) ( 3p 7 + 28p 5 + 178p 3 + 28p + - ) arctanp - 13440p 4 In 


p 2 + 1 


(Cll) 


It was mentioned in Sec JTT] that in order to obtain the ’’pure” component V’i' 2 ? one needs 
to select and then to get rid of admixture of the HH, Y&. In other words, we need to 
orthogonalize to each of Y^. First of all, according to expansion (193|) and Eqs. (fTTf - (fT2|l 
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we need to calculate the unnormalized HH expansion coefficients 

•7r/2 


Fu = n 2 Nl 


da t[ v \p){ sin a) 


1+2 


Y 4 i(a, 6 )Pi(cos 6 ) sin 9d9+ 


da r/ p ' ) (l/p)(sina 


/ 71-/2 


, 1+2 


> 4 ;(a, 9)Pi(cos9 ) sin 9 d9 


(C12) 


where the factor N 41 is defined by Eq.Q. The direct integration with the use of solutions 
flmP-flCHl) yields: 

(tt-2) 


Fin = 


81007T 2 

F 41 = 0, 

>-2) 


[5tt( 15 In 2 - 16) + 247 - 300G], 


F42 — 


[5tt( 840 In 2 + 109) + 4592 - 16800G], 


(C13a) 

(C13b) 

(C13c) 


113400tt 2 

where G is the Catalan’s constant. Thus, to obtain the ’’pure” subcomponent one 

needs to substitute by r 0 = — J- 4 0 ( 4 cos 2 a — 1 ), and r^ by r 2 = r^ — F 42 in 

expansion (l93lh 

For l > 3, the particular solution of Eq. (lC5l) obtained by Eq.(j45]), can be presented in 
the form: 


Ti P) (P ) = 


1—2 


~ 2)(p 2 + 1) 

1357r2 3 0+ 1 )A;p 2i+1 


1—3 


P 


21+2 


dnP n + PY1 P^P 211 ~ 4o (P 2 + !) 6 arctal1 P 4 


(3) 2 n 


n =0 


n =0 


n =0 


(C14) 

where the coefficients A 1 and are presented in Table UT1 whereas p Y and p^ are presented 
in Table ITTT1 in the form of lists. The number of coefficients is limited by l < 10 . 

It is clear that, all the terms of expansion (|93li with l > 3 are automatically orthogonal 
to the Y 41 . Therefore, in this case we shall use relation (lOTTh which becomes: 


7 21,1 — 


(/ + !)! 


b(a)(sina) 2 Z+ 2 da: = 


yJYY(l + 3/2) J o 


n(p)- 


P 


21+2 


2 2 0+ 2 )(7 + 1 )! r 1 

(i + 3 / 2 ) A (i + p2)2<+3 


-dp, 


(C15) 

where F n j are the unnormalized HH expansion coefficients for the considered subcomponent 

4 2 ?(m) = XXwm)- (° 16 ) 

nl 

The prime indicates that n = 4 is omitted from the summation for ’’pure” component. 

The general solution of Eq. (lC5ji reads 


Tl(p ) = r/ P) (p) + AiUuip) + A 2V4 i(p), 


(C17) 
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where the solutions and v<u of the homogeneous equation associated with Eq. (lC5D are 
dehned by Eqs. Phjl . and the constants A\, A 2 are currently undetermined. It is clear that in 
order if}ff 1 satisfies the finiteness condition P2, function (sin a) l Ti(p) for each / must satisfy 
this condition, according to expansion (jU3|) . It is easy to verify that (sin a) l u±i(p ) is divergent 
at p = 0 (a —> 0), whereas both functions (sin a) l v^i(p) and (sina)V^(p) are finite at this 
point. Hence, one should put A\ = 0. Thus, Eq. (lC15|) can be rewritten in the form: 

Tv = +iT(/+3/2) lVl(l) + ' (C18) 


where 


v,(i) = 


v 2 (l) = / n 4 i{p)- 


P 


Tl P \p ) 

21+2 


P 


21+2 


-dp, 


(1 + p 2 ) 2l+3 

n 2~ 2 ^T (l + 3/2) 


-dp = 


(C19) 


(C20) 


(1 + p 2 ) 2i + 3 r r(z/2 + 3)r(z/2) 

Note that for each value of l > 3, the integral (1C 1911 can be calculated in the exact explicit 
form. 

To derive the closed expression for T n ,i according to definition (1C16|) . first, let us obtain 
the HH expansion for the rhs of Eq. flCip . 




(C21) 


nl 


where according to Eq. (IT2]) for unnormalized HH, 


n ( ‘) = K, I h^’(c,e)Y nl (a,e)dn, 


(2d) 


(C22) 


and the normalization coefficient N ni is dehned by Eq.(|8]). Substitution of representation 
(IC2h - (IC3h into (1C22h yields: 


-i/(4) _ 
H n,l ~ 


KA* ~ 2 ) 


3(2 1 - l)(2l+l)2 l ~ 1 J o 
5 


sm a 


, 21+2 


COS 


a 


sm 


a 


Ci^lfcosct) x 


2 J w n/2—l 


FiAa) - ( 1 


Fi,i(a) - (21 - l)Fi_i(a) 


da. (C23) 


(21 — 3) sin a ' ' ' \ sma^ 

The orthogonality of the Legendre polynomials was used to derive the last formula. Inserting 
expansion (1C16(1 into the lhs of Eq. (lCl|) . one obtains: 


(A 2 - 32) %,i(n - 4)(n + 8 )Y nl (a, 9). 


(C24) 


nl 
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According to Eq. (lCl|) . the right-hand sides of Eq. (1C24|) and Eq. flC21|) can be equated, what 
yields: 


T n ,l = 


n 


(4) 

n,l 


(n — 4)(n + 8) 

Thus, for the requested case of n = 21 (l > 2), one obtains, using (1C25I) and (1C23I) : 

(n -2)2 l ~ 1 (l+l)(l\) 2 


T 21 J — 


3vr 2 (2/ — !)(/ —2)(/ + 4)(2/ + l)! 


w), 


(C25) 


(C26) 


where 


)2Z+2 

cos 

(a 

\2 

^ + sin 

Fi,s(a) - 

0 

2 

sin a 


x 


Fi,i( a ) — (2/ — 


(21 — 3) sin a 

Integral (1C27[1 can be taken in the explicit form. The result is 

2 1+1 


da. 


(C27) 


V 3 (l) = 


(l + 3) (21 - 3) (21 + 3) (21 + 5) 
30 26 


x 


+ 13 — 4/ [47 - 21(21(1 + 3) - 9)] + 2 ;+1 x 


l+l 1+2 


3 1 


(/ + 1) [4Z(Z(4Z + 3) - 17) + 45] Bi ( l + 4 J + 8/ [1(1(41(1 + 4) + 3) - 56) - 62] Bi ( l + ^ 


where B z (a,b) is the incomplete beta function. Equating the right-hand sides of Eq. (IC18|) 
and Eq. (lC26p . one obtains the required coefficient in the form 

(71 - 2 ) 2 ~ 3 6 + 2 ) 


A 2 = A 2 (l) — 


;V 3 (l)-V l (l) 


(C29) 


V 2 (l) |_37t(2Z — 1)(Z — 2)(l + 4) 

It may be verified that for odd values of /, the coefficients A 2 (l) equal zero. For even values 
of l. formula (IC29|) yields 

(rr — 9 ^ 

(C30) 


am = -++2a(i), 


7 


where 


,, N 5515vr- 11648 ., . 191095vr - 396032 . 66779345 tt - 137592832 

•A(4) = , “4.(6) = , “4(8) = - 


2948400 


x4(10) = 


1378377000 
59227659tt - 121716736 

25014766104900 ' 


4085509428000 


(C31) 


(C28) 
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Appendix D 


The IFRR (T75|) for subcomponent reads 

(A 2 - 21) = C. 


(Dl) 


where the rhs h^ is dehned by Eq. (l88]h To find a suitable solution in the form of expansion 
), let us present the rhs of Eg. (ID ID in the form 


h (2c) = 

n 3,0 ~ 


4£ 


3 sin a 


Pi( cosP)(sina) i h/(a), 


(D2) 


/=o 


where, using expansion (1A7D for v = 1, one obtains 

, , i 2 2_z ^ (l 11 1 , 3. 2 

= TTTo;— TV~- - 2 ^M 7 _ 7’ o + 7’ + o’ sm a 

3(2/ —l)sma \2 4 2 4 2 


(D3) 


Inserting expansions (1D2D and (1981) into the IFFR (1D1D . equating then the factors for the 
Legendre polynomials P;(cosP) of the same order, and turning to the variable p = tan(a/2), 
one obtains the inhomogeneous differential equation 


(l + p 2 ) <f>i(p) + 2p 1 [l + p 2 + /(I — p 4 )] 4>i{p) + (3 — 21) (7 + 2 l)4>i(p) — —hi(p), (D4) 

which presents Eq. (l44|i with k = 3, and with changing h(p) for hi(p) = lp(o;), and g(p) for 
0 z (p). Simplification of the rhs of Eq. (JD4[) yields for 0<P<1: 


h{p) = 


)1-Z ( „2 


(p 2 + 1)' + ^ [(1 - 2Z)p 2 + 2/ + 3] 


(D5) 


3(2/ - l)(2/ + 3)p 

Simplifying the solutions (146|) of the homogeneous equation associated with Eq. (1D4|) . one 
obtains for k — 3: 


u 3 i{p) = 


(P 2 + 1) 


i-i 


(21 + 3) (2/ + 5)^ 4 + 2(2/ + 5) ^ 2 + 1 


vM = {p 2 + l) 


p 2i+i [(2/ — 3) (2/ — l) r ' 2/ — 1 

(21 — 3)(21 — 1) 4 2(2/-3) 


z-# 


:P + 


P 2 + l 


(D6) 


(D7) 


(21 + 3) (2/ + 5) ' 2/+ 3 

The particular solution ) of Eq. (lD4D can be obtained by formula (j45jh but with changing 
definite integration for indefinite integration (that is putting the antiderivatives for the lower 
limits equal zero). The result is presented by Eqs. (1101 j) - (ll04j) . It is worth noting that the 
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rypergeometric function included into Eq. (11041) can be presented in the form (see, 7.3.1.135 

22 ]): 


2 Fi (l, l + 1; l + 2; —p 2 ) — 


(l + 1 ) 


-P' 


2V+ 1 


In (1 + p 2 ) + ^ 


2\ k 


k =1 


(V) 

k 


(D8) 


It can be easily shown that the homogeneous solution u 3 i(p) is singular, whereas v 3 ;(p) and 
the particular solution 0^(p) are regular at the point p = 0 (a = 0) for any / > 0. The 
conclusion is: the physical solution must be sought in the form (1991) . The problem is to find 
the coefficient q. To solve the problem we shall apply the coupling equation (16TT) . which for 
this case becomes: 


021,1 — 


(/ + 1)!2 2 ( ;+2 ) 7 1 

0Fr (/ + |) ,/o 


fi\p) + C t V 3 l(p) 


p 


21+2 


J (P 2 + 1) 


21+3 


dp , 


where O n j present the unnormalized HH expansion coefficients for subcomponent 




a,0) - ^TlVjWo.O). 


(D9) 


(DIO) 


nl 


For this case (odd value of k — 3) the rhs of Eq. fIDlOp presents the physical solution we 
are looking for (see, the end of SecJTTl). To derive the closed expression for O n j, first, let us 
obtain the unnormalized HH expansion 

h^\ a ,e) = Y, H n}Yn,(a,e) 


(Dll) 


nl 


for the rhs of Eq. dDlj) . where by definition 


H^=Nl] h^(a,e)Y nl (a,e)dn. 

Inserting expansion (1D2|1 - (1D3D into ( II ) 121 ) , one obtains 


,(2c) 


(D12) 


" nl " 7 - / 2 Fi ( - — sin 2 a ) (sino;) 2 ' +i cos a)da, 

n ’ 1 3(2/ — 1)(2Z +1) J 0 V 2 4 2 4 2 J 7 

(D13) 

where the orthogonality of the Legendre polynomials was used. Substitution of expansion 
into the lhs of Eq. (ID 11) yields: 

(A 2 - 21) V4 2 0 c) = ]T O nil (n -3)(n + 7 )Y nl (a, 6). 


N 2 nl n 2 2 3 -* 


/ 1 / 


\2£+l /^i{l~\~ 1) 


(D14) 


nl 


According to Eq. dDl|h the right-hand sides of Eq. dDllD and Eq. dD14D can be equated, what 
yields: 


O n , = 


n 


(3) 

n,l 


(n — 3)(n + 7)' 


(D15) 
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Thus, inserting (ID 13[) into Eq. (ID15|) for n = 21 (l > 0), one obtains: 

2 l+4 l\(l + 1)! 


021,1 ~ 


t/2 


„ , l 11 1 . 3. 2 . 

2 -f i [ - - -, - + Z + sm a ] (sin a 


2 4 2 4 

2 2 ~ l l\(l + 1)! 


3tt(2/ - 3)(2/ - 1)(2Z + 7)(2 1 + 1)! 
21+1 da = 


x 


'Z 1 Z 1 , „ , 3 , 3 „ 

3-^2 - — Tj ~ + Tj Z + 1; Z + -, Z + -; 1 


3(2/-3)(2/-l)(2/ + 7)T 2 (/ + 3/2) V 2 4 7 2 ' 4'" ' 7 ' ' 2 7 " ' 2 

Equating the right-hand sides of Eqs. (ID9|) and (lD16|h we find the required coefficient, 

M 1 (l)-M 2 (l) 


(D16) 


Cl = 


Ms(l) 


(D17) 


where 

Mi(Z) = 


>—3Z—2 




3(2 1 - 3)(21 - 1)(2Z + 7)T(Z + 3/2) 


„ .'21 — 1 21 + 1 , , 3 , 3 , 

3-^2 ( -7-,-7-, l + 1; l + 77, l + 775 1 > 


•Ad 3 (Z) = / V3i(p) 


p 


2i+2 


4 ’ 4 

y-l- + 


(P 2 + 1) 


2£+3 


dp = 


2-‘-s(2Z + l) 
(2/ + 3)(2Z + 7)' 


(D18) 

(D19) 


According to representation (ri01l) - (U04[) . function A4 2 (l) can be presented in the form 


MM = I +\p) 


p 


21+2 


(P 2 + 1) 


2i+3 


dp = 


3(21 — 3)(2Z — 1)(2Z + 3)(2Z + 5) 


2A4 21 (Z) + 


2M 22 (l) + A4 23 (Z) 
21 + 1 


where 


M 21 (l) = - 


6 


13 — 4Z 2 8Z 3 + 28Z 2 — 21 — 79 /, 7 9 , 

-— +-—- 2^1 ( Z + 2, Z + -; Z + 3; —1 


2 l+ i 


1 + 2 


(D20) 


(D21) 


M 22 (l) = 


(2Z + l)(2Z + 5) r (21 + 3) 


(2Z + 3)(2Z + 7) 2 ( 2 l+ 2 

2^(1 + 1 )! 


vr(2Z + 7) 


4(4Z 2 + 24Z + 23) 
(2Z + l)(2Z + 5) _ 


. + 2a (Z + 1) 2 Fi ( — — Z; - 

r(z + |) 1 ; V 2 2’2 


M 23 (l) = (-1) 


/+1 




k =1 


(- 1 )* 
k 


Ikl 


(D22) 

(D23) 


In the last function we have introduced the following notations: 

2“*-i 


V, = 


181 x 2 l+ * - 525 In 2 - 1046 + Z 


2 /+ ^~ (l + 3) (11 + 21(1 + 3)) 


(2Z + 3)(2Z + 5)(2Z + 7) 2 

-8(195 + 4Z(31 + Z(Z + 9))) - 2(985 + 4Z(303 + 2Z(83 + Z(21 + 21)))) In 2 


(D24) 
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Qu = (2 ' 2( ^ 2 ' 1) 11 2 *i (k + 1, (+ |; fc + 2; -l) + (2i + 5) x 

/ 9 \ 2/ + 3 / 9 ^ 

Fl ( /c + 2 ’ /+ 2 ;/C + 3;_1 J + 2(fc + 3) 2jPi ( fc + 3 ^+2 ;fc + 4; - 1 y 


21 -3 


/c + 2 


To derive formulas (1D23[) ~ (1D25[) . the representation (ID8D was used. 


. (D25) 
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TABLE I: The simplest subcomponents of the solutions, and the corresponding rhs of the FRR. 


k 

P 

3 

rt k,p 


1 

0 

0 

-2/C 

C/2 

3 

1 

1 


1 

Uy> 

1 CO 

^co 

3 

0 

0 

Ei+ ^r 1 

c I" 1-2 E , (S-2) t i 
“C 24 1 72 “C 

3 

0 

la 

5(w-2)t 

3tt s 

5(tt-2) .3 

72tt s 

3 

0 

2 a 

-K2C+C- 1 ) 

(1 - £0 

4 

1 

1 

^ [6(1 - 2 E) + (12 E - 5)C 2 ] 

|gg[3(32E-15)-8(12E-5)C 2 ] 

3 

0 

16 

^[^ + 2(1-3E)] 

[4 E — 1 + (E — 1) sin a] 

3 

0 

3 


— ^ (2 + 5 sin a) ? 

4 

1 

2 b 

2, "S- U> (3 4sin 2 a) 

^~ 540 ^~ 14 ^ csc a [ acos (3ce) — g sin(3a)] 

4 

1 

3 

[2 + tan(a/2) + cot(a/2)] sin a cos 6 

— (4 + 5 sin a) sin a cos 9 


TABLE II: Numerical coefficients /+[[ (columns 2-6) and A i included into Ea. (IC14l) . 


l\n 

0 

1 

2 

3 

4 

A i 

3 

-522720 

-1051985 

255552 

-90435 

227872 

495 

4 

-3953664 

-8398510 

3727360 

-419475 

3335168 

2002 

5 

-30351360 

-97464675 

44154880 

-6969375 

39739392 

6552 

6 

-2882764800 

-9137296644 

5778063360 

-189629055 

5234966528 

235620 

7 

-270052392960 

-1196759060145 

696562155520 

-65054175615 

635252473856 

7759752 

8 

-2853364039680 

-12357527011830 

9051962343424 

-67609913085 

8305992794112 

27387360 

9 

-5615652962304 

-31944952238815 

21219319480320 

-1350136582575 

19579564720128 

17341632 

10 

-132056403148800 

-733168078648720 

580357875302400 

1595198130669 

538208290996224 

127481640 
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TABLE III: Numerical coefficients fij^ and /t® included into Ea. (jC14p . 


l\!i 

(2) 

MnZ 

(3) 

Vnl 

3 

{-72765, -412335, -960498, -2625678} 

{1} 

4 

{1273965,5945170,9597203,3639900, -14747117} 

{1,-1} 

5 

{-8513505, -34999965, -45846801, -10135125,10135125, -146271983} 

(l,-f,l| 

6 

{1456717185,5473725180,6250641012,983757060, -777756650, 

983757060, -13174180236} 

{l,-fl,fl,-l} 

7 

{-101948591745, -358127104335, -370912265724, -46034728740,31617756170, 

-31617756170,46034728740, -1585784814084} 

fi 28 378 28 
l 1 ’ 13’ 143’ 13’ 

i} 

8 

{1804845357315,6016151191050,5790931172103,602937350136, -374632125450, 

328821494780, -374632125450,602937350136, -16476813837561} 

r-i 7 42 42 

t 1 ’ 3’ 13’ 13’ 

1,-1} 

9 

{-2508790869585, -8018292387105, -7290251115147, -663280843995,381499104294, 

-306553004730,306553004730, -381499104294,663280843995, -39837938484597} 

fi 42 63 924 

l ’ 17’ 17’ 221’ 

63 42 i \ 

17’ 17’ 

10 

{95748321954285,295643941472880,256791268415682,21004496987760, 

-11353909326045,8538680495904, -7845295711140,8538680495904, 

-11353909326045,21004496987760, -931653655078718} 

r i 49 1323 1617 
l 1 ’ 19’ 323 ’ 323 ’ 

1617 1323 49 11 
323 ’ 323 ’19’ J 
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